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ANALYSIS METHOD AND REAL TIME MEDICAL OR COGNITIVE 
MONITORING DEVICE USING THE ANALYSIS OF THE CEREBRAL 
ELECTROMAGNETIC ACTIVITY ON AN INDIVIDUAL, APPLICATION 
OF THIS METHOD TO CHARACTERISE AND DIFFERENTIATE 
BETWEEN PHYSIOLOGICAL AND PATHOLOGICAL STATES 

DESCRIPTION 

Technical domain 

This invention relates to an analysis method and a 
device for real time medical or cognitive monitoring 
using the analysis of the cerebral electromagnetic 
activity of an individual, with applications of this 
5 method for differentiation and characterisation of 
physiological and pathological states, particularly for 
real time processing of anticipation of epilepsy 
- seizures , 

10 State of prior art 

Role of cerebral interactions in man: dynamic mapping 

Any cerebral action is the result of a cooperation 
between several spatially distributed neuron networks. 

15 At the present time, and despite recent progress, the 
principal cerebral imaging techniques, namely EEG 
(electroencephalography) , MEG (magnetoencephalography) , 
FMRI (functional magnetic resonance imaging) and PET 
(positron emission tomography) , only provide a map of 

2 0 cerebral activations without directly reporting 
interactions between these activations. 

Characterisation of these functional networks requires: 
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- identification of the cerebral regions involved; 

- understanding of interaction mechanisms between 
them ; and 

- precise quantification of these interactions. 

5 It is impossible to observe operation of these 

networks simply from a map of cerebral activities. How 
would it be possible to decide which zones among all 
active zones at any one time participate in a 
particular network? A simple observation that two 

10 zones are active at the same time is not enough to 
conclude that they are engaged in the same pathological 
or cognitive process. To elucidate these exchange 
mechanisms, it .is necessary to have explicit 
assumptions about the nature of these links and to have 

15 technical means of observing them. 

All approaches known in prior art are based on the 
idea that the existence of a coupling between two zones 
must result in a correlation between their activities. 
However, the activity of a group of neurons, for 

20 example a cortical column, can be characterised by two 
measurement types : 

- time coding with the rate of neuron discharges 
per second; or 

coding by synchronisation of oscillatory 
25 activities of cerebral zones involved in the same 
network . 

The purpose of the invention is to propose a 
method for dynamic mapping of the brain starting from 
such a phase synchronism measurement, starting from the 
30 assumption according to which dynamic links between 
neuron groups occur by synchronisation of oscillatory 
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activities in some frequency bands between 0 and 
2000 Hz. 

Example application to pathology: anticipation of 

5 epilepsy seizures 

In the remainder of the description, the method 
according to the invention is described with reference 
to medical monitoring of a patient for anticipation of 
epilepsy seizures in real time, simply as an example 

10 and to make the description clearer. Obviously, it 
would be possible to apply the method according to the 
invention to other fields and to characterise and then 
recognise physiological or pathological states by 
changing the database . 

15 Epilepsy is one of the most frequent neurological 

complaints in children and adults (1% of the 
population) , and is the result of neuron disorder 
expressed by recurrent paroxysmal or paroxystic 
discharges of the cerebral cortex. The clinical result 

2 0 of epilepsy is the sudden occurrence of the symptoms of 
a seizure. This sudden emergence is difficult to 
interpret as a response to a triggering external 
factor, which is not present in most situations, except 
for the rare reflex epilepsies. The transition between 

25 an " intercrit ical " state and the critical state (the 
seizure) is one of the primordial phenomenologies of 
epilepsy and this intermittence appears to be an 
unpredictable product of an internal self -organization 
of the brain. 



4 



No traditional method (linear analysis) has yet 
been able to significantly anticipate this state of 
seizure . 

Two recent publications, references [1] and [2] at 
5 the end of the description, describe a process that 
under some conditions can make it possible to 
anticipa:te seizures lasting several minutes using new 
strategies output from the theory of dynamic systems. 
Non- linear dynamic methods are derived from known 

10 mathematics under the generic term "Chaos theory" . 
They show how there can be precise laws or determinisms 
behind an apparently random elect roencephalographic 
signal. The possibility of anticipating the occurrence 
of seizures opens up very broad prospects. 

15 Anticipation of seizures would be useful for the 

large number of patients with a drug resistant 
epilepsy, which accounts for about 25% of all epileptic 
patients. The intermittent and unexpected occurrence 
of seizures is a recognised factor of mortality (by car 

20 accidents or sudden death) and morbidity (particularly 
brain and face traumatisms) . Furthermore, patients 
consider the unpredictable nature of seizures to be one 
of the most important factors affecting their poor life 
quality. Limitations related to the risk of seizures 

25 are social (isolation due to the fear of a seizure in 
public) , professional (risk activities that epileptics 
cannot perform) and personal (car driving prohibited) . 
The possibility of anticipating the occurrence of 
seizures would be a means of helping epileptics to live 

30 better with their seizures. In addition to the 
possibility of alerting the patient about the imminence 
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of a seizure, anticipation strategies should also be 
helpful for making complementary examinations for the 
purposes of a pre- surgical checkup of their epilepsy. 
Furthermore, such strategies should open up. the 
5 opportunity for therapeutic actions before the seizure 
has time to emerge. 

Document reference [1] mentioned above is based on 
quantification of the similarity between a "normal" 
reference dynamic and the emerging epileptic dynamic. 

10 This similarity index is calculated independently for 
each recorded channel. The space and time components 
of cerebral changes can be obtained by aligning changes 
in the statistical deviation of the indexes of each 
channel. Nevertheless, this method only indirectly 

15 takes account of the spatial aspect of modifications to 
electroencephalography variations (intracranial and 
surface) that we use to predict the .imminence of a 
seizure. However, strong evidence suggests that this 
period reflects the transition from a disordered state 

2 0 to a more ordered (or less complex) state that would 
very likely correspond to changes in the. 
synchronisation of several neuron assemblies 
distributed in different cerebral structures. The 
ideal method of increasing the efficiency of 

25 anticipation would be to have additional information 
about the space-time behaviours of epileptogenic 
variations such as subtle variations in activity or 
interdependences between distant cerebral regions. 
Many observations suggest that it is often very 

30 difficult to determine a single zone of epileptic 
malfunctions as being the origin of a seizure. In 
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particular, recent work has very much suggested the 
importance of imagining partial epilepsies as being a 
manifestation of a networked structure. It is quite 
probable that propagation of the discharge facilitates 
5 a number of connections between multiple neuron 
assemblies, thus more efficiently modifying local and 
remote neuron connections. Thus, the organization of 
epileptic malfunctions can no-, longer be envisaged 
statically or as a local malfunction (circumscribed 
10 epileptic focus) , but rather satisfies a complex 
dynamic space-time model that spatially involves neuron 
networks connected by abnormally facilitated 
connections involving some synchronised behaviours in 
time. 

15 Linear techniques (cross-correlations in the time 

domain or consistencies in the frequency domain) have 
frequently been used in the past to study the spatial 
organization of the epileptogenic network and to 
characterise interactions between the epileptogenic 

20 network and the remainder of the brain. In the case of 
an analysis between macroscopic signals, these methods 
are often limited by assumptions made in them related 
to the stationarity pf signals and the linear nature of 
interactions. This is particularly problematic because 

25 non-linear behaviour is significantly increased in an 
epileptic period. 

The purpose of this invention is to overcome these 
limitations . 
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Presentation of the invention 

The invention relates to a method for analysing 
synchronisations of the electroencephalography of an 
individual using a set of sensors starting from the 
5 cerebral electromagnetic analysis of the patient, 
characterised in that it comprises the following steps 
consisting of: 

- a step to create a database comprising: 

• a phase for acquisition and digitisation of 
10 electrophysiological signals output from these sensors, 

• a phase to calculate the degree of 
synchronisation existing between all pairs of sensors 
recorded in an assembly protocol, in frequency bands 
between 0 and 2000 Hz, to build up this database of 

.15 classes each characterising a reference state; 

- a step for statistical validation of a period 
analysed in real time, which assigns this period to a 
class in the database, 

a step to detect a specific period with a 
20 determined degree of synchronisation. 

Advantageously, ^ the said method includes an 
analysis associated with at least one type of 
electrophysiological signals among electrocardiograms, 
electrooculograms, electrodermograms, breathing 

25 signals. 

Advantageously, a PLS method is used during the 
statistical validation step, which estimates the phase 
difference between oscillations of signals from two 
electrodes. The statistical level of the PLS 

3 0 synchronisation between two signals is evaluated using 
the circular variance of the phase difference between 



8 



the signals or using the normalised Shannon entropy of 
the phase difference between the signals. 

The method according to the invention can be used 
to characterise and differentiate between phyisiological 
5 or pathological states, for example for anticipation of 
epilepsy seizures . 

The method according to the invention can be used 
in other application fields, such as: 

- sleep: differentiation between different sleep 
10 stages; 

anaesthesia: characterisation of stages of 
falling asleep under anaesthesia with automatic control 
of regulation of the injected substances- 
depression: with electrophysiological 
15 monitoring of a depressive patient and characterisation 
of his traits or states and consequently adjustment of 
his treatment; 

schizophrenia: with electrophysiological 

monitoring of a patient and quantification of his 
20 traits or states for diagnosis assistance and 
therapeutic purposes; 

- diagnosis assistance for neurological diseases 
such as Parkinson's and Alzheimer's diseases; 

- characterisation of cognitive states (levels of 
25 vigilance and attention, perception and conscious and 

unconscious recognition of visual, auditory, 
somes thetic and emotional stimulations (fear, joy, 
etc . ) . 

The invention also relates to a real time medical 
30 or cognitive monitoring device starting from the 
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cerebral electromagnetic analysis of an individual, 
characterised in that it comprises: 

means of acquiring and digitising 
electrophysiological signals output from sensors; 
5 - means of calculating the synchronisation between 

all pairs of sensors recorded in an assembly process, 
in frequency bands between 0 and 2 000 Hz, to build up a 
database of classes each characterising a reference 
state; 

10 - means of statistical validation of a period 

analysed in real time to assign this period to a class 
in the database; 

means of detecting a cognitive period or a 
specific pathological period; 

15 - means of sending an alert signal if applicable. 

Advantageously, the device according to the 
invention is a standalone, lightweight device that the 
patient can carry himself or herself. The device 
according to the invention may be miniaturised so that 

20 it can be implanted subcutaneous ly like a stimulator, 
so that patients can be completely autonomous. 

Brief description of the figures 

- Figure 1 illustrates the different steps in the 
25 method according to the invention, starting from the 

electroencephalogram (EEG) analysis . 

Figure 2 illustrates the steps in the method 
according to the invention more precisely. 

- Figure 3 illustrates the device according to the 
30 invention. 
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Figure 4 illustrates explanatory chronograms 
showing the treatment of electroencephalograms using 
the method according to the invention. 

- Figure 5 illustrates an example embodiment of 
5 the method for real time treatment of epilepsy- 
according to the invention. 

Detailed presentation of particular embodiments 

Neuroelectric activity in a restricted frequency 

10 band is characterised by its energy and its phase, such 
that demonstration of a relation between two groups of 
neurons within a given frequency band requires the 
demonstration of a significant correlation between 
variations in the energy or phase of each. The most 

15 frequently used method at the moment makes simultaneous 
use of the energy and the phase. It consists of 
calculating consistency between signals using the 
"Magnitude Squared Coherence" (MSG) index. 

This MSG index is a global measurement in which it 

20 is difficult to separate the influence of the phase 
from the influence of energy. A correlation between 
variations of the phases of the two signals may be 
sufficient to demonstrate coupling between two neuron 
groups (document reference [3]) . 

25 The method according to the invention enables a 

measurement of the synchronism, using the phase only, 
namely the "Phase Locking Statistics" (PLS) method. 
For a given latency, this method estimates the phase 
difference between oscillations of signals from two 

30 electrodes. If this phase difference is relatively 
constant during the analysed period, a high PLS index 
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will be obtained which is a sign that the two 
electrodes are significantly synchronous. 

This PLS method is sufficiently precise to detect 
periods with synchronism and is therefore suitable to 
5 describe a series of transient synchronisms like those 
presumed to occur in cognitive processing or to 
characterise more sustained synchronisms like those 
presumed to characterise pathological states. 

This method can be used to measure the degree of 

10 synchronisation between the activities of various 
cerebral regions. Synchronism between two neuron 

groups is defined within a given frequency band as 
being a significant correlation between variations in 
their phase with time; this is called "phase-locking" . 

15 Nevertheless, considering volume effects (the activity- 
of a single neuron population can thus be picked up by 
two relatively remote electrodes) and neuron background 
noise, the only way to detect synchronism between two 
regions is to use a statistical approach. The 

20 statistical validity of the measurements is then tested 
by the construction of bi -varied surrogate data. 

The invention can thus use a statistical 
estimating method based on the use of surrogate data, 
which the PLS method uses to apply itself to 

25 non- stationary neuroelectric signals (unlike the PSC 
method), as is the case for most biological signals. 

Therefore the method according to the invention is 
a real time medical or cognitive monitoring process 
based on an analysis of the cerebral electromagnetic 

30 activity of an individual associated with analysis of 
any other electrophysiological signals 
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(electrocardiograms electrooculograms, 

electrodermograms, breathing signals), particularly to 
detect specific cognitive or pathological periods, for 
example an epilepsy seizure under preparation, and to 
5 provide an alert signal in all cases necessary to enable 
prevention or a therapeutic action. 

As illustrated in Figure 1, the process according 
to the invention includes the. following steps: 

- an electrophysiological signals acquisition and 
10 digitisation step 10: in general, a headset placed on 

an individual's scalp and provided with 27 to 128 
electrodes depending on the problem to be solved, can 
be used to record the individual's cerebral activity 
with a fairly good spatial resolution. Some additional 
15 signals may be acquired at the same time (eye movement 
signal, heart activity signal, etc.); 

- a step 11 to calculate synchronisation between 
all pairs of signals and in. several frequency bands to 
build up a database (step 12) of reference states 

20 depending on the selected problematics (pathological, 
sleep, wake, vigilance, etc.); 

- a step 13 for statistical validation of the 
period analysed in real time, to classify this period 
starting from the database. This validation is based 

25 on a non-parametric multidimensional discrimination 
method; 

a step 14 to detect specific cognitive or 
pathological periods ; 

a step 15 for sending an alert signal, if 
3 0 applicable. 
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More precisely, the following steps can be 
performed in sequence starting from a database with k 
classes, as illustrated in Figure 2: 

- calculate the partition of the variables space 

5 by Bayesian probabilities S^^^f ; 

- sort a time window x, for example lasting for 10 
seconds, in the base of k classes; 

- detect a class with alert if applicable. 

As illustrated in Figure 3, the independent real 
10 time medical or cognitive monitoring device using an 
electromagnetic cerebral analysis of an individual 
includes circuits (amplifier 20, analogue-digital 
converter 21, buffer 22) for acquisition of signals 
representing the electrical activity of the brain, a 
15 processor 23 being used for acquisition and processing 
of these signals and an alert circuit 24 for the 
patient or for his environment, for example a light. 

Expected results in epilepsy and clinical 

20 implications 

It has been observed that some pairs of electrodes 
in the periphery of the epileptogenic zone 
systematically show a significant modification of their 
degree of synchronism before a seizure, particularly 

25 for example in the alpha (8-12 Hz) , beta (15-30 Hz) and 
gamma (30-70 Hz) frequency bands. Interestingly, this 
synchronisations have recently been given considerable 
attention for their possible role in large scale 
integration phenomena during cognition. These results 

30 thus suggest that neuron populations subjacent to the 
epileptogenic zone modify their relations with larger 
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scale variations before the seizure. These changes in 
synchronisations may cause "dynamic isolation" of the 
focus and could recurrently provide a neuron population 
that can easily be recruited by epileptic processes. 
5 New analyses techniques for synchronisations of 

electroencephalography used in the method according to 
the invention can be used to very precisely quantify 
the pre-critical cerebral activity. This possibility 
of anticipating the occurrence of seizures opens up 
10 very broad medical prospects: 

- in fundamental research, by characterisation of 
neurobiological modifications that occur during this 

. precritical phase; 

- clinically, by the possibility of warning the 
15 patient, and attempting to stop the initiating seizure 

by a therapeutic action. 

In particular, electric neurostimulation recently 
appeared as a promising therapeutic solution for other 
pathologies, particularly for Parkinson's disease. In 

20 this context, ^ the mechanical . destruction of a 
predefined cerebral region has been replaced by the 
principle of a protective treatment by electric 
stimulations to reinforce or inhibit a neuron activity. 
In this respect, the possibility of anticipating 

25 seizures using the method according to the invention is 
decisive since it provides a solution to the question 
of "when to stimulate?". These stimulations may be 
applied when a pre-ictal state is detected and are 
aimed at destabilizing epileptogenic processes before 

30 they become irreversible at the time of the seizure. 
This is the approach that is preferred in the long term 
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for patients being investigated by intracerebral 
electrodes . 

Furthermore, the analysis technique used for 
synchronisations of electroencephalography can enable 
5 developments of "cognitive" actions. Some patients 
describe their ability to interrupt their initiating 
seizure by carrying out specific cognitive activities 
or by motor activities. These phenomena are probably 
based on destabilisation of the epileptic process by 

10 the appearance of new electrical activities within the 
cerebral cortex. Inventors have also demonstrated 
modulation of an epileptic activity by cognitive 
synchronisations . 

Other actions could also be applied, for example 

15 pharmacological action consisting of administration of 
fast active antiepileptic drug (like benzodiazepines) . 

These possibilities of alert and actions available 
if seizures can be anticipated, necessarily require 
"real time" anticipation, in other words the results of 

20 mathematical calculations need to be obtained 
instantaneously in the method according to the 
invention, not later on. 

The ability to anticipate seizures can also 
improve examinations performed during the pre -surgical 

25 check-up of partial pharmacoresistant epilepsies. In 
particular, the pre-critical cerebral . scintigraph 
(SPECT-ictal) is facilitated if the team is alerted: 
injection of the radioactive tracer immediately before 
or at the beginning of the seizure improves the 

30 definition of the location of the epileptogenic centre. 
Hospitalisation times can then be considerably 
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shortened and imaging system occupancy times can be 
optimised. The possibility of anticipating the 

occurrence of epilepsy seizures due to in depth and 
surface electroencephalography opens up a fairly broad 
5 range of prospects in social and clinical applications. 

Mathematical procedure used for calculating phase 
synchronism between two signals 

The instantaneous phase of a signal may be 
10 calculated either using an analytic signal, a concept 
introduced by Gabor in 1946 and recently applied on 
experimental data, or by convolution with a complex 
specific wavelet (Lachaux et al. Human Brain Mapping, 
1999) . 

15 For an arbitrary signal s(t), the analytic signal 

^ is a complex function dependent on time defined as> 

follows : 

^(t) = s(t) + j?(t) = A(t)e^^f^^ (1) 
where the function J (t) is the Hilbert transfoirm 
20 of s (t) : 

s {t)=-Py^^dT (2) 

P.V. states that the integral is calculated in 
terms of the Cauchy principal value. The instantaneous 
amplitude A(t) and the instantaneous phase (j)(t) of the 

25 signal s(t) are only defined by equation (1) . As can 
be seen in equation (2) , J (t) is considered to be the 
convolution product of s(t) and I/tc. This means that 
the Hilbert transform is equivalent to a filter for 
which the amplitude response is unitary and the phase 

30 response is offset by n/2 for all frequencies. 
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Although this transform can be applied in theory to 
signals with a wide frequency band, the phase concept 
is not very clear in this case. In practice, only 
narrow band signals obtained by filtering are, used. 
5 Consequently, filtering is always done in a specific 
frequency band. Several frequency bands may be 

selected. The same frequency band is used for two 
signals in case of 1:1 synchronism. Different 
frequency bands are used for the study of mm 
10 synchronisms. The statistical level of PLS 

synchronisation between two signals is evaluated using 
one of the following two indexes: 

- the circular variance of the phase difference 
(A<|)) between signals; or 
15 - the Shannon normalised entropy of this phase 

difference (A<|)) . 

The circular variance is such that: 
M 

Shannon's normalised entropy is such that: 
20 Y=(Hn,ax - H)/H 

max 

where the entropy H is defined by: 

H = EPkl^Pk 
k=l 

where N is the number of classes, Hniax=ln(N) is the 
maximum entropy, and pk is the relative frequency of 
25 the phase difference in the k^^ class. The optimum 
number of classes is N=exp[0 . 626 + 0 . 41n (M-1) ] where M is 
the number of elements (phase difference) to be sorted. 
With this normalisation, the values of y for 95% of 
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"surrogates" (replacement values) are included between 
0 (uniform distribution and no synchronisation) and 1 
(perfect synchronisation) . This calculation is made 
for all pairs of sensors recorded in the assembly 
5 protocol. The number of distinct pairs for the 27 
electrodes in a standard assembly is 26x25/2=325, and 
the number for 128 electrodes is 8001. 

Thus, as illustrated in Figure 4, we have the 
following in sequence: 
10 - passband filtering of two signals obtained with 

two electrodes 3 0 and 31 (f ± 1 Hz) ; 

- Hilbert's transform of these signals; 

- evaluation of the statistical level of the PLS 
synchronisation using two indexes: 

15 • entropy of Acj) (phase difference between <j)l and 

(|>2) ; 

• circular variance of A<t>, 

The second step consists of setting up the 
database of calibrated states of the patient as a 

20 function of the objective to be achieved. 

The third step is a discrimination step for 
decision-making purposes. Considering a recording 
period equal to 10 seconds in pathology, but sometimes 
much shorter for the discrimination of cognitive states 

25 for which quantification by the synchronisation method 
is known, an attempt is made to assign it to a class 
characterising one cerebral state among several. This 
is a sort problem, which assumes that a set of classes 
has been defined in advance. The main difficulty is 

30 the dimension of the variables space. Quantifying 
synchronisation between all pairs of sensors in six 
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frequency bands requires a variables space with 
dimension p=1950 (325x6) for a setup with 27 
electrodes. The expression for the a posteriori 
probability of the analysed time window x belonging to 
5 the k different groups of cerebral states is as follows 
(Bayes's theorem) : 

K 

P(Gr/x)=P(Gr) .P{x/Gr)/ Z P (Gj ) , P (x/Gj ) where r = 1, 

where P(Gr) is the a priori probability of 
10 belonging to a class and in practice is estimated by 
the frequency of elements of Gr in the total sample. 
The different values P (x/Gr) are estimated by 
probability densities. For each new analysed period x 
to be sorted into one of the k groups, a search is made 
15 for the q closest neighbours of each of the k groups, 
thus defining the average radius rk of the hypersphere 
HS(rk,x) containing the average of the q neighbours of 
x and the volume of the corresponding hypersphere in 
the space R^. Thus, the probability density P(x/Gr) 
20 can be estimated as follows: 

P(x/Gr) = — ^ 
n^Aj- 

and X is assigned to the group je[l,k] if 
:P(Gj/x)=maxr{p(Gr/x) ; r=l,2,..., kL 

25 Example use of the invention 

The device according to the invention was applied 
to electrical intracranian cerebral recordings of 
patients being sxibjected to surgical treatment of their 
temporal epilepsy and it has thus been demonstrated 
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that seizures can be anticipated by several minutes and 
there is a deterministic "road towards seizure" 
phenomenon, as illustrated in Figure 4. The device 
according to the invention enables the seizure in the 
5 10-20 Hz frequency band to be anticipated almost 
20 minutes before the seizure, and is characterised by 
reduced synchronisation. The power spectrum, a 

classical signal processing procedure, does not give 
such clear modifications. 
10 The method according to the invention is 

particularly suitable for clinical situations and makes 
it possible to extend the results to surface 
electroencephalography, due to its low sensitivity to 
recording artefacts • 
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